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Abstract 

We present an antenna shower formalism including contributions from initial-state partons and correspond- 
ing backwards evolution. We give a set of phase-space maps and antenna functions for massless partons 
which define a complete shower formalism suitable for computing observables with hadronic initial states. 
We focus on the initial-state components: initial-initial and initial-final antenna configurations. The formal- 
ism includes comprehensive possibilities for uncertainty estimates. We report on some preliminary results 
obtained with an implementation in the Vincia antenna-shower framework. 
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1. Introduction 

Parton-shower algorithms offer a universal and 
fully exclusive perturbative resummation frame- 
work for high-energy processes. In the context of 
Monte Carlo event generators Q], they also pro- 
vide the perturbative input for hadronization mod- 
els. As such, they are complementary to more 
inclusive techniques, such as fixed-order calcula- 
tions (limited to small numbers of hard and well- 
separated partons) and more inclusive resumma- 
tion approaches (limited to a fixed set of observ- 
ables). 

Sjostrand derived the first consistent parton- 
shower algorithm for so-called "backwards 
evolution" of initial-state partons a quarter-century 
ago. The central point is that an initial-state par- 
ton defined at a high factorization scale, Qp , can 
be evolved "backwards", towards earlier times, to 
find the parton from which it originated at some 
low scale, Qq ~ 1 GeV. During this evolution, 
which is governed by the Altarelli-Parisi splitting 
kernels Q supplemented by PDF ratios (a point 
which is crucial to the backwards-evolution for- 
malism), initial-state radiation is emitted, which in 
turn gives rise to its own final-state radiation, and 
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the character of the evolving parton changes, mi- 
grating towards successively higher x values and 
towards the more valence-dominated flavor content 
at low Q. 

As an alternative to Altarelli-Parisi evolution, 
Gustafson and Pettersson proposed a final-state al- 
gorithm based on QCD dipoles H, which has 
been implemented in Ariadne Q. There, how- 
ever, initial-state radiation does not rely on back- 
wards evolution. Instead, it is treated essentially as 
final-state radiation off dipoles stretched between 
the hard process and the beam remnants, and thus 
depends on the non-perturbative makeup of the 
remnants. Winter and Krauss took a first step to- 
wards combining the dipole formalism with back- 
wards evolution (and thus also eliminating the de- 
pendence on the remnants) in ref. [6]. Our con- 
struction differs in the antenna functions, evolution 
variables, and recoil strategy. In particular, it dif- 
fers in the treatment of collinear singularities in 
initial-final antennae. We have checked that our an- 
tennae properly reproduce all QCD singularities. 

A complementary approach which merges the 
Lund dipole language with that of fixed-order an- 
tenna factorization Q HJ HI [TO), is that of Vin- 
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Figure 1 : Illustration of initial-initial and initial-final branch- 
ings: AB — > a jb and AK — > a jk, respectively. For the II case, 
the recoil of the hard system is illustrated by the change in 
orientation of the three outgoing lines representing the origi- 
nal final-state system. 



cia 031 Q21 EES. (Note: we henceforth use the 
term "antenna" rather than "dipole" to avoid am- 
biguities of historical origins, see e.g., ref. [14]). 
So far, however, the Vincia formalism has been ap- 
plied only to final-state showers. In this paper, we 
present all the ingredients necessary to construct a 
consistent initial-state shower based on QCD an- 
tennae. A further important ingredient is compre- 
hensive possibilities for uncertainty estimates, in 
line with the framework for automated theory un- 
certainties proposed in ref. lPT5l . 

2. Antennae and Antenna Showers 

Throughout this Letter, we use the following no- 
tation convention: capital letters for pre-branching 
(parent) partons and lower-case letters for post- 
branching (daughter) ones. Also, we use the first 
letters of the alphabet, a,b,c, for incoming par- 
tons and letters starting from i, j, k,... for outgoing 
ones. Fig. [T] illustrates these choices for the two 
basic types of configurations we consider. We will 
also mark incoming particles with a minus sign 
in front in antenna functions. We adopt the con- 
vention that particle energies are always positive, 
whether the particle is in the initial or the final state. 
As a result, sij = (ki + kj) 2 is always positive. 

The key building block for parton showers is the 
Sudakov factor, which represents the non-emission 



probability between two values of the evolution 
scale, see CD HH for reviews. In the context of an 
antenna shower, the Sudakov factor for the branch- 
ing of one antenna is 
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In this equation, dO ant represents the antenna 
phase-space factorization, which provides an ex- 
act Lorentz-invariant mapping from 2 to 3 on-shell 
partons, that conserves global energy and momen- 
tum. Specific forms appropriate to initial-final and 
initial-initial antenna configurations are defined in 
sections [3] and [4] respectively. 

The evolution variable Q 2 is a function of the 
phase-space point and must vanish in the unre- 
solved limits [17]. The general formalism permits 
us to study different evolution variables iTTTl [T5ll , 
though in this Letter we will restrict ourselves to 
a transverse-momentum type variable, defined in 
section [5] As in all parton showers, the descrip- 
tion is expected to be accurate only in the strongly- 
ordered limit for the Q 2 of successive emissions. 

The dressed or colored antenna function a c is de- 
fined a£] 



a c - AnasiQ )Ca , 



(3) 



where C is a color factor (we recall that we use nor- 
malization conventions such that gluon and quark 
emission antennae have C = Ca and C = 2Cp, re- 
spectively, and gluon-splitting ones have C = 1), 
and a is a color-ordered antenna function, which 
embodies the factorization of QCD matrix ele- 
ments in all single-unresolved soft and collinear 
limits. We don't take the functions a to be fixed; 
instead we use different antenna functions with the 
same singular limits as one estimate of the shower 
uncertainty. 

1 Note that in [ 15 1 the normalization was a c = a s l{An)Ca 



We use so-called global antenna functions @ 
(called sub-antenna functions with uniquely identi- 
fied radiators in ref. [9]) which are active over all of 
phase space. A backwards-evolution shower based 
on sector antennas in analogy to refs. (T8] [13l is 
left for future work. Some, but not all, antennae 
needed for initial-state radiation can be chosen to 
be the crossings of their final-final counterparts. 
An incoming particle is necessarily a hard radiator 
in an antenna. Therefore, a gluon emission antenna 
function with an incoming gluon has to reproduce 
the AP splitting function on its own, e. g. 

1 
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whereas if both gluons are in the final state, the 
collinear singularity is reproduced by the sum of 
two antenna functions 
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where the first antenna function is singular for i be- 
coming soft, the second for j. 

In pure final-state showers, the x values of the 
incoming partons are not modified by the phase- 
space factorization, hence the PDF ratios in eq. ([2]) 
drop out, yielding the ordinary form of the final- 
final Sudakov form factor QT1[T5). 

For initial-final antennae, only one of the PDF 
x values changes, and a Sudakov factor very 
similar to that of conventional AP showers re- 
sults, with a single PDF ratio in the kernel, 
fa(xa, Q 2 )//a(xa, Q 2 )- Unlike conventional show- 
ers, however, we must also consider the back- 
wards evolution of two initial-state partons simul- 
taneously, generally requiring two separate parton- 
density factors in initial-initial antennae. 

The consideration of initial-initial and initial- 
final antennae gives rise to one more subtlety. The 
basic antenna functions are color-ordered, so that 
in a final-final gluon-emission antenna, for exam- 
ple, the emitted gluon is color adjacent to both 
other (hard) daughter partons. That is, it is the 
middle parton of the color trio which is emitted. 
The leading-color approximation inherent in par- 
ton showers along with the symmetry of final- state 



phase space allows us only antennae with this order- 
ing. When considering initial-state antennae, how- 
ever, the emitted parton need not be color- adjacent 
to both other daughter partons; the middle parton, 
adjacent to both, may end up in the initial instead 
of the final state. We will call antennae in which 
the middle parton is emitted into the final state, 
'emission' antennae; and those in which the mid- 
dle parton ends up in the initial state, 'conversion' 
antennae. 

For those antennae in which the type (spin) of the 
initial-state partons does not change after branch- 
ing, we can redistribute collinear singularities to 
neighboring antennae so as to replace 'conversion' 
antennae by 'emission' antennae. For those an- 
tennae in which the type of the initial-state par- 
tons changes during branching — in which a quark 
backwards-evolves into a gluon or vice versa — we 
cannot avoid a consideration of both types of an- 
tenna function and non-emission probability. 

3. Initial-Final Configurations 

The pre- and post-branching partons for initial- 
final configurations are labeled by AK — > a jk, with 
the other incoming parton, B, acting as a passive 
spectator, see the illustration in fig. [T] 

In general, the incoming momentum after 
branching will no longer be parallel to the beam di- 
rection. We could boost it back to the beam direc- 
tion; this will transfer some of the transverse mo- 
mentum generated in the emission to the rest of the 
event. This is the antenna analog of the recoil con- 
sidered in ref. [fT9l . In the present Letter, we will 
instead restrict the branching so that the incoming 
momentum remains parallel to the beam axis after 
branching. 

With this restriction, the phase-space factoriza- 
tion reads [20], 
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with xaIx ci = SAK/(sAK + Sjk) and where the initial- 
final antenna phase space is 
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with the boundaries < Sjt < sak(1 - *a)/*a, < 
s a j < SAK + Sjk- We have suppressed the integration 
over the third coordinate of the initial-final phase 
space on which the emission probability does not 
depend. 

The gluon-emission antennae can be chosen as 

• , x 1 (ZSakSAK s jk Saj\ ... 
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where the factor 1/2 originates from the fact that 
the gluon is part of two antennae. 

The antenna governing the backwards-evolution 
of a gluon into a quark is 

7 7 
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For the reverse process of a sea quark 
backwards-evolving into a gluon, we use 
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with a color connection j — a — k at variance with 
the other antennae. 
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where it is apparent that the antennae with an 
incoming quark are crossings of their final-state 
counterparts whereas the ones with incoming glu- 
ons have additional terms compared to their final- 
state counterparts to ensure the collinear singular- 
ity a || j is taken into account properly. 

The antenna for the splitting of a final- state 
gluon into a quark-antiquark pair is chosen as 
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4. Initial-Initial Configurations 

For initial-initial antennae, we label the pre- and 
post-branching partons by AB — > a jb, see fig. [T] 
In the initial-initial case, we must necessarily have 
transverse momentum generated, which must be 
absorbed by the rest of the event. There are two 
ways of proceeding. One can allow the incom- 
ing partons to be shifted away from the beam di- 
rection after branching, and then boost back to a 
frame in which they are again parallel to the beam 
direction. Alternatively, one can fix the incoming 
partons to be parallel to the beam direction, and 
balance the new transverse momentum by boosting 
the rest of the event appropriately. In both cases, 
there is a freedom in how the longitudinal part of 
the emission momentum is absorbed into the ini- 
tial state. This corresponds to a freedom in relat- 
ing the post-branching momentum fractions x a> b to 
the pre-branching momentum fractions xa,b- I n tne 
first case, this freedom is parametrized by the recoil 
or reconstruction function r in combination with 
the Lorentz transformation boosting back to the lab 
frame. In the second case, it is parametrized by the 
functional form of x a j,- 

It turns out that these two approaches are equiva- 
lent, unlike the initial-final case. We define our re- 
coil strategy in terms of x a ,b here. The phase-space 



4 



factorization reads lETTl . 
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with the initial-initial antenna phase space 
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where we have suppressed the integration over the 
angle <p parametrizing rotations around the beam. 
The pre- and post-branching momenta are related 
by a Lorentz transform: 
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The phase space boundary depends directly on 
the definition of the post-branching momentum 
fractions, which is not fixed completely by the re- 
quirements of x a XbSAB = XAXBS a b and the behavior 
in the soft and collinear limits. 

For gluon emission, we use 
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whereas for conversion, we keep one incoming mo- 
mentum fixed, i.e. x A lx a - sabI s a b, Xb = xb, giv- 
ing the phase space boundaries s a j + sjb ^ *ab(1 _ 
xa)/xa- This corresponds to the use of a one-sided 
factorization, which is possible because only one 
collinear limit is singular, whereas the behavior of 
the phase space factorization in the other collinear 
limit is not constrained. 

We use the emission antennae 
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For a quark backwards-evolving into a gluon, we 
use 

1 / 2SjbS A B ^ S a b 
3 aj\Sab ~ Saj) S a j 

(22) 

The antenna for a gluon backwards-evolving 
into a quark is the crossing of the initial-final coun- 
terpart: 
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5. Implementation and Preliminary Results 

In the antenna shower, as in a conventional 
shower, we start the evolution at high Q 2 , and gen- 
erate a series of branchings at successively lower 
Q 2 , stopping when we reach a shower cut-off, typ- 
ically around 1 GeV. Each branching is generated 
according to the non-emission probability ([T]), and 
in this Letter we shall restrict ourselves to strict 
strong ordering, postponing a discussion of smooth 
ordering [ 15 ] and/or power showers |[22ll to a sub- 
sequent study. 

In order to generate a branching, we must in- 
vert the function specified by the integral Q. This 
is in general a difficult task even if the integral 
is doable analytically, because the result involves 
dilogarithms. In some cases, the boundaries even 
make it unreasonable to perform the integral ana- 
lytically. A direct inversion would in either case 
be quite slow. Instead, we proceed as follows. 
We pick a simple function — a trial antenna func- 
tion a 1 " 3 ' and trial ratios of parton-density func- 
tions — which overestimates the integrand, 



and veto the excess emissions generated accord- 
ing to the non-emission probability computed us- 
ing the trial function. The trial function is cho- 
sen to capture the leading logarithmic singularities 
of the antenna function, and to allow the phase- 
space integral to be factorized into a product of 
one-dimensional integrals. Where possible, it is 
also chosen to produce an analytically invertible in- 
tegral. In the final-final case, the latter requirement 
can always be satisfied; in initial-initial and initial- 
final cases, it can be satisfied for most trial an- 
tennae. In the exceptional cases, we employ a two- 
stage veto. In these cases, the first-level trial func- 
tion still serves to simplify the inversion of the non- 
emission probability by ensuring the factorization 
of the integral into a product of one-dimensional 
integrals. The veto probability is given by, 
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evaluated at the post-branching kinematic point. In 
this equation, 
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is the ratio of partem densities that appears in 
eq. Q. 

When approaching a heavy flavor threshold, /? pt jf 
diverges. As is standard in backwards-evolution 
codes, we absorb the leading divergent behavior of 
Rpdi into the trial integration to maintain a reason- 
able efficiency Il23l . Analogous issues may arise 
at large x and low Q 2 in light-quark parton densi- 
ties due to numerical instabilities. We defer their 
treatment to future work. 

To define a concrete shower algorithm based on 
the above antennae and phase-space factorizations, 
we have chosen to use two different evolution vari- 
ables, depending on the type of antenna. For gluon 
emission, we use a transverse momentum, 
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For final-state branchings, Q 2 ± is equal to 2p 2 ±A , the 
evolution measure used in Ariadne [5] (note: pre- 



vious Vincia publications used 4p 2 ±A ). The max- 
imal value of Q 2 ^ in the final-state case is Sijk/2. 
For conversion and gluon splitting antennae, we in- 
stead use the virtuality of the only potentially sin- 
gular propagator as the evolution variable. As in 
the final-state shower ifTTl [121 [131, other choices 
are possible within the Vincia formalism. We de- 
fer an exploration of more general possibilities to 
future work. 

We now turn to a few basic tests of each com- 
ponent of the shower algorithm. In all cases, we 
consider pp collisions at 8 TeV CM energy, use 
the MSTW 2008 LO PDF set [24], with a one- 
loop running as , normalized to as(niz) = 0.13939. 
In all calculations performed here, we turn off 
hadronization and primordial kj, as well as the un- 
derlying event, both in the Vincia calculation and 
in the Pythia 8 [25 ] calculations to which we com- 
pare. While the evolution variables in Vincia and 
Pythia are different, we have tried to match the 
shower cut-offs in calculations with the latter to 
that we use in Vincia. This includes accounting 
for the difference in normalization between lim- 
iting definitions of transverse momentum. Note 
that while Vincia uses a zero-mass variable flavor 
number scheme, Pythia 8 uses the physical quark 
masses everywhere. For the observables we dis- 
cuss here, the effect is negligible. 

In fig. [2] we show the pj spectrum of the Z 
boson in Drell-Yan production, which is sensitive 
to radiation in initial-initial configurations. The 
main figure pane shows the peak of the distribution, 
while the inset shows the high-pr tail. The figure 
shows Vincia curves computed using two different 
antenna functions: the default antenna given ear- 
lier, and an enhanced antenna function, with a finite 
term — 5/ Sjju — added. It also shows the result 
obtained with Pythia 8. 

The overall shape of the three curves is similar: 
small values at small pj, rising to a peak and then 
declining again with a rough power-law fall in the 
asymptotic region, and a "knee" around pj ~ 90 
GeV due to the requirement of strong ordering 
which we have imposed here. The difference be- 
tween the two Vincia predictions illustrates the un- 



0.03 



■jjj 0.02 
as 



0.01 



20 40 60 80 100 120 



°0 




■ Pythia 

■ Vincia (default) 
Vincia (enh. antennae) 



10 

P 



15 



20 



Figure 2: The Drell-Yan p T spectrum. The dashed red curve 
shows the value computed using Vincia with default antennae 
functions, while the dotted green curve shows the Vincia pre- 
dicted with an enhanced antenna function. The solid blue 
curve gives the Pythia 8 prediction. The inset shows the high- 
p T tail. 



certainty due to the shower function and in particu- 
lar higher-order terms in the shower. The differ- 
ence shown here is illustrative only; a more ex- 
tensive exploration of possible antenna variations 
would be required before taking the spread as a 
quantitative estimate of the uncertainty. We may 
nonetheless observe that the Pythia 8 reference 
calculation differs from the Vincia one (with de- 
fault antenna) by roughly the same amount in the 
peak region as does the enhanced Vincia predic- 
tion. This illustrates a tradeoff between a more ac- 
tive recoil strategy (Pythia) and a more active radi- 
ation pattern (enhanced Vincia), which will be in- 
teresting to study more closely. At large pj, all 
three curves are close to each other; the transverse 
momentum here is dominated by the recoil against 
hard lone-gluon emission. This region would be 
described well by fixed-order calculations. 

For initial-final configurations, coherence is par- 
ticularly important, and can lead to sizable asym- 
metries (see, e.g., Il26l ). An illustration of the effect 
is given in fig. [3] which shows qq — > qq scatter- 
ing with two different color-flow assignments: for- 
ward (left) and backward (right). In both cases, 
the starting scale of the shower evolution would 
be pj, the transverse-momentum scale character- 
izing the hard scattering. Coherence, however, im- 




Figure 3: Different color flows and corresponding emission 
patterns in qq — > qq scattering. The straight (black) lines are 
quarks with arrows denoting the direction of motion in the ini- 
tial or final states, and the curved (colored) lines indicating the 
color flow. The beam axis is horizontal, and the vertical axis 
is transverse to the beam. The initial-state momenta would be 
reversed in a Feynman diagram, so that the gluon emissions 
symbolically indicated by curly lines would be inside the cor- 
responding color antennas. Forward flow is shown on the left, 
and backward flow on the right. 
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Figure 4: Angular distribution of the first gluon emission in 
qq — > qq scattering at 45°, for the two different color flows. 
The light (red) histogram shows the emission density for the 
forward flow, and the dark (blue) histogram shows the emis- 
sion density for the backward flow. 



plies that radiation should be directed primarily in- 
side the color antenna, so that in the forward flow 
it would be directed towards large rapidity, and 
strongly suppressed at right angles to the beam di- 
rection. In the backward flow, conversely, radiation 
at right angles to the beam should be unsuppressed. 
The two radiation patterns are illustrated schemat- 
ically by the gluons in fig. [3] The intrinsic coher- 
ence of the antenna formalism accounts for this ef- 
fect automatically. That Vincia reproduces this fea- 
ture is demonstrated in fig. |4j which shows the an- 
gular distribution of the first emitted gluon for the 
forward and backward color flows, respectively, for 
a scattering angle of 45° and pj = lOOGeV. The 
distributions clearly show that the backward color 
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Figure 5: The dijet decorrelation angle. The histograms are 
normalized to unity separately. 

flow allows for much more radiation at 90° than the 
forward one. The pj spectrum of the radiation (not 
shown) is also harder. The next step will be to inter- 
face the hadronization and underlying-event mod- 
els in Pythia, and compare to experimental studies, 
such as the one by CDF ETl (we note that an up- 
date of that study, correcting it to the hadron level, 
would be highly useful to the MC community). 

Finally, to demonstrate the combination of all 
shower components acting together, we show the 
dijet decorrelation angle (the azimuthal angle be- 
tween the two leading jets), A<pjj, in fig. [5] Using 
FastJet l28l . we consider anti-^T jets with radius 
parameter R = 0.4. We demand the two leading jets 
to have transverse momentum above 100 GeV and 
to be at rapidities \y\ < 2.8. Note that Pythia 8 pro- 
duces about 40% more jet events which pass these 
cuts than Vincia both with or without enhanced an- 
tenna functions, partly due to its more active re- 
coil strategy, which allows the original dijet system 
to build up transverse momentum successively dur- 
ing the shower cascade. This difference is not visi- 
ble in fig. [5J as the distributions are all normalized 
to unity. The two Vincia distributions are broadly 
similar to the Pythia 8 distribution, and very sim- 
ilar in the two-jet region (Acpjj ~ n) where the 
parton-shower approximation should be reliable. 
The differences are substantial in the region below 



A<pjj < 3n/4, where hard real emission is impor- 
tant. In this region, fixed-order calculations should 
be reliable, but unmatched parton showers will not 
be. Nonetheless, the difference between the Vincia 
calculation with default antenna strength and the 
Pythia 8 calculation is similar to that between the 
two Vincia calculations, suggesting again that the 
variation provides a good qualitative assessment 
of the uncertainty. We defer a comparison of this 
distribution with fixed-order calculations to future 
studies. 



6. Conclusion and Outlook 

In this Letter, we have presented the outline of a 
formalism for an antenna shower for hadron col- 
lisions, along with results from an initial imple- 
mentation as a plug-in to Pythia 8. The formal- 
ism requires the introduction of new antenna?, cor- 
responding to one or both parents being initial- 
state partons. These should be further subdivided 
into the new categories of emission and conver- 
sion antennae based on their color flows. The 
formalism also requires factorizations suitable for 
phase spaces involving initial-state partons, and in- 
troduces ratios of parton densities into the non- 
emission probabilities governing the shower evo- 
lution. 

We have chosen to implement the shower as 
a plug-in to the Pythia 8 program, which takes 
advantage of the latter's flexible framework and 
makes use of its utilities, structures, and overall 
management of the branching process. In this ap- 
proach, it replaces Pythia's shower with an antenna 
shower. For practical and efficiency reasons, we 
uniformly adopt a trial-and-veto algorithm for gen- 
erating branchings. The trial functions used in the 
implementation will be described elsewhere. 

We expect to implement further optimizations of 
the branching step in future work. The leading- 
order matching approach described in ref. lfl~5l 
should carry over to the initial-state showering de- 
scribed here, and will be an important next step for 
the development of Vincia. 
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